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Preface 


This  technical  report  contains  part  of  the  work  completed  under  the  in-house  project  work 
unit  2301  SI  10,  project  title  “Volume  and  Surface  Plasmas”.  Some  of  the  earlier  work 
was  published  in  technical  report  #  WRDC-  TR-90-2027.  The  project  engineer  for  this 
task  is  Dr.  Biswa  N.  Ganguly.  The  research  for  this  work  was  performed  in  the  period  01 
February,  1989  through  01  September,  2004. 


I.  Introduction 


This  report  describes  work  performed  to  obtain  (1)  scaling  rules  for  the  generation 
of  low  pressure  nonequilibrium  plasma  for  applications  such  as  anisotropic 
plasma  etching  of  wide  band  gap  semiconductors,  surface  passivation  of 
semiconductors  devices  for  improved  operating  lifetime  under  high  radiation  flux 
environment;  (2)  energy  efficient  atom  flux  generator  through  multiquantum 
vibrational  energy  transfer  in  gas  mixture  plasma;  (3)  remediation  of  dust 
particles  in  rf  plasma  processing  equipment;  (4)  self-assembled  carbon  nano¬ 
particle  production  in  both  reactive  and  non  reactive  plasmas  through  sputtering; 
(5)  In-situ  plasma  and  sheath  electric  field  measurements  in  both  volume  and 
surface  dominated  plasmas  by  Rydberg  State  Stark  Spectroscopy;  (6)  absolute 
atom  density  measurements  by  two-photon  allowed  laser  induced  in  dc,  rf,  and 
microwave  excited  plasmas;  (7)  actinometric  measurements  of  reactive  atom 
density  in  rf  etching  plasmas;  (8)  strong  double  layer  formation  by  weak  shock 
propagation  through  nonequilibrium  plasmas;  (9)  shock  mitigation  by  plasma- 
shock  interactions. 

The  main  objectives  of  this  project  has  been  to  develop  methods 
for  quantitative  measurements  of  plasma  properties  which  may  permit  to  develop 
a  scaling  law  for  a  specific  application  such  as  high  flux  atom  source  generator  or 
surface  damage  free  ion  assisted  etching  of  semiconductor  materials  with  large 
cohesive  energy.  Similarly,  for  quantification  of  plasma-  aerodynamic 
interactions,  an  electrically  driven  shock  tube  instrument  was  developed  which 
pennit  simultaneous  measurements  of  both  gas  dynamic  and  plasma  properties 
modulated  by  the  plasma-shock  interactions.  This  approach  has  allowed  us  to 
quantify  the  relative  efficiency  of  shock  mitigation  by  plasmas  through  direct  and 
indirect  gas  heating.  The  details  of  some  of  the  recent  research  results  are 
described  in  sections  II,  III,  IV,  V  and  VI.  Also,  other  previous  published  research 
results  that  were  published  in  the  referred  journal  publications  and  as  invited 
paper  presentations  are  listed  in  the  Appendix  section.  This  in-house  task  has 
resulted  into  one  issued  patent  one  patent  pending.  Also,  more  than  twenty  journal 
publication  and  two  book  chapters  have  been  written  based  on  the  research  results 
obtained  under  this  in-  house  task.  Four  Ph.  D.  theses  were  also  sponsored. 
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II.  Two  color  Mie  scattering  images  of  uncorrelated  and  correlated 
trapped  dust  in  low  frequency  helium  discharges. 


We  report  measurements  of  self  assembling  of  dust  particles  by  their  size 
distribution  using  simultaneous  near  real  time  two  color  Mie  scattering  in  a  2  Torr  pure 
helium,  15  to  20  KHz  voltage  excitation  frequency,  discharge  operated  with  graphite 
electrodes.  A  3  cm  i.d.  and  1  cm  wide  cylindrical  graphite  is  used  as  outer  electrode  and 
a  1  cm  diameter  graphite  is  used  as  center  electrode  which  are  mounted  in  a  5  cm 
diameter  glass  tube  (shown  in  figure  1).  In  this  discharge  configuration,  the  dust  particles 
are  transported  by  ion  drag  force  to  an  equilibrium  position  just  outside  the  interelectrode 
volume  and  they  are  trapped  by  the  radial  electric  field.  This  condition  leads  to  a 
minimum  perturbation  of  the  discharge  electrical  characteristics  as  the  dust  particle 
density  grows.  In  pure  helium  discharge  with  graphite  electrodes,  the  particle  growth  is 
primarily  by  ion  accretion  of  sputtered  carbon  cluster  ions.  Ion  accretion  usually  favors 
particle  growth  with  a  narrow  size  distribution1.  This  characteristic  of  the  discharge 
permit  us  to  follow  the  time  evolution  of  trapped  particles  by  their  size  distributions  and 
also  the  transition  of  the  particles  from  uncorrelated  phase  to  a  phase  with  long  range 
correlation.  A  krypton  ion  laser  output  of  647.  lnm  and  an  argon  ion  laser  output  514.5 
nm  are  combined  using  a  broadband  half  silvered  mirror  so  that  each  beam  is  propagating 
parallel  to  each  other  but  spatially  displaced  to  permit  chopping  by  the  alternate  blades 
using  a  mechanical  chopper.  A  beam  expander  combination  lens  is  used  for  both  beams 
to  illuminate  a  spatially  overlapping  wide  discharge  volume.  A  color  video  camera  is 
used  at  90  degree  angle  to  detect  the  Mie  scattered  light  from  both  laser  beams  at 
alternate  video  frame.  A  video  frame  grabber  board  is  used  for  data  transfer  to  a 
computer.  Since  the  Mie  scattering  efficiency  ~  scales  as  (d/A)  ,  where  d  is  the  particle 
diameter  and  X  is  the  scattered  light  wavelength,  at  a  fixed  scattering  angle,  the  smaller 
particles  will  preferentially  scatter  514.5  nm  whereas  the  larger  particles  will  scatter 
647.1  nm.  This  arrangement  permits  imaging  of  self-assembled  dust  particle  growth  in 
both  uncorrelated  and  correlated  phases.  Fig.  1  and  Fig.  2  (left)  show  the  images  of  the 
trapped  dust  particles  at  the  edge  of  the  plasma  after  run  time  of  several  minutes.  Note 
that  larger  particles  (red  scattered  light)  are  assembled  closer  to  the  electrodes  compared 
to  the  smaller  particles.  Also,  the  smaller  particles  are  trapped  and  distributed  over  a 
larger  volume  compared  to  the  heavier  particles.  The  spatial  variation  of  scattering 
intensities  of  both  647.1  nm  and  514.5  nm  show  the  dust  densities  are  concentrated  on  or 
near  the  axis  of  symmetry  with  respect  to  the  electrodes.  Under  this  discharge  operating 
condition,  there  is  no  measurable  amount  of  dust  particle  trapped  in  the  interelectrode 
volume  due  to  the  relatively  strong  radial  electric  field  variation  which  permits  transport 
of  the  negatively  charged  particles  to  the  boundary  walls. 

After  run  time  of  tens  of  minutes,  the  Mie  scattering  intensities  of  both  laser  beam 
undergo  a  transformation  from  diffuse  distributed  light  scattering  to  crisp  objects  with 
very  well  delineated  boundaries,  corresponding  to  the  formation  of  highly  correlated  dust 
ensemble.  Fig.  1  and  Fig.  2  (right),  along  with  the  blow  up  insets,  show  the  Mie  scattering 
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images  of  both  scattered  lights.  These  dust  ensembles  are  located  near  the  axis  and 
correspond  to  the  previous  locations  of  highest  dust  densities,  in  the  uncorrelated  phase, 
which  implies  the  presence  of  a  very  small  radial  electric  field  gradient.  Also,  the  dust 
particles  of  same  size  are  self-  organized  to  form  correlated  ensembles. 

Besides  the  characteristics  of  highly  delineated  boundaries,  these  dust  ensembles  exhibit 
very  different  dynamic  behavior  compared  to  the  uncorrelated  dust  particles.  The 
transport  properties  of  these  correlated  dusts  differ  significantly  from  the  uncorrelated 
phase  during  the  discharge  switch  off.  The  uncorrelated  dust  particles  undergo  Coulomb 
explosion  with  the  dispersal  of  the  particles  in  all  directions  including  against  gravity. 
Whereas,  the  dust  ensembles  in  the  correlated  phase  retain  their  crisp  boundaries  for 
several  video  frames  (>  0. 1  sec)  and  their  motion  is  determined  primarily  by  the  balance 
between  the  Stokes  force  and  the  gravitational  force.  This  persistence  suggests  that  in 
such  dust  ensembles  the  strongly  correlated  3  ion-  ion,  and  electron-ion  motions  lead  to  a 
several  orders  of  magnitude  reduction  in  the  ion-ion  and  electron-ion  recombination  rates 
on  dust  particles  compared  to  the  ambipolar  or  to  free  diffusion  conditions. 

References: 

1.  P.  Haaland,  A.  Garscadden,  and  B.  Ganguly,  Appl.  Phys.  Lett.  69,  904  (1996). 

2.  H.C.  van  de  Hulst,  “Light  scattering  by  small  particles”(Dover,  New  York,  1981),  p 
144. 

3.  V.N.  Tsytovich,  U.de  Angeles,  R.  Bingham,  and  D.  Resendes,  Phys.  Plasmas  4,  3882 
(1997). 
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Fig  2. 

Figure  Captions: 

Fig.  1 :  Mie  scattering  image  from  uncorrelated  (left)  and  correlated  (right)  dust  using 
5 14.5  mn  laser  beam. 

Fig. 2:  Mie  scattering  image  from  uncorrelated  (left)  and  correlated  (right)  dust  using 
647. 1  mn  laser  beam. 
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ITT.  The  effect  of  displacement  current  on  fast-pulsed  dielectric  barrier 

discharges. 


1 .  Introduction 

Dielectric  barrier  discharges  (DBD)  are  in  wide  use  for  ozone  generation,  toxic  gas 
remediation  and  surface  treatment  of  materials.  A  conventional  DBD  uses  high  AC 
voltage  at  power  line  or  higher  frequency.  The  configurations  may  include  a  parallel  or 
cylindrical  electrode  configuration  with  one  or  both  electrodes  covered  with  dielectric 
material.  The  breakdown  of  the  gas  under  these  conditions  typically  consists  of  a  series  of 
microdischarges  at  the  positive  and  or  negative  peak  of  the  applied  voltage.  The  exact 
nature  of  these  microdischarges  is  determined  by  the  gas  pressure,  type  of  gas  and  the 
repetition  frequency  ’  '  of  the  applied  voltage.  With  certain  gases  such  as  He  or  N2  and 
with  the  applied  peak  voltage  only  slightly  above  the  gas  breakdown  voltage,  the 
discharges  may  be  limited  to  low  currents,  and  this  type  of  discharge  may  be  in  the  form 
of  an  atmospheric  glow  discharge.  The  breakdown  voltage  at  the  operating  frequency  of 
the  applied  AC  voltage  determines  the  reduced  electric  field  E/n  (where  E  is  the  electric 
field,  n  the  neutral  gas  density).  Due  to  the  buildup  of  surface  charges  on  the  dielectric 
barrier,  the  E/n  decreases  during  the  discharge3  and  may  also  decrease  with  consecutive 
micro-discharges.  Plasma-chemical  reaction  rates  (for  processes  such  as  excitation, 
dissociation)  increase  rapidly  with  increasing  E/n  and  it  is  therefore  desirable  to  operate 
the  DBD  at  a  high  E/n  for  high  plasma-chemical  efficiency.  It  has  been  recognized  4’ 5 
many  years  ago  that  the  high  E/n  could  be  achieved  by  operating  the  DBD  with  pulsed 
excitation.  For  example  Lowther  4  used  a  10  kHz,  0.8  ps  6kV  pulsed  DBD  and  claimed 
an  ozone  production  efficiency  of  more  than  three  times  of  that  of  a  high  voltage  corona 
discharge  with  oxygen  feed  gas  and  almost  six  times  that  of  dry  air  feed.  Rosocha  5  used 
a  parallel  plate  Blumlein  pulser  that  also  incorporated  the  discharge  cell  to  investigate 
ozone  production  with  pulses  of  up  to  75  kV  amplitude  and  risetimes  as  fast  as  5  ns. 
Measurements  and  calculations  of  the  effect  of  E/n  in  pulsed  corona  or  pulsed  DBD’s 
show  different  results,  in  some  cases  because  of  different  gas  species.  Some  experimental 
results  show  an  increase  of  chemical  reactivity  with  E/n  4’6,7'8  and  while  others  show  no 
effect  of  increasing  E/n  9’10,11. 

The  breakdown  voltage  using  a  pulsed  voltage  with  short  pulse  width  is  normally 
higher  than  when  using  a  slowly  rising  DC  voltage,  since  the  probability  of  an  ionization 
event  due  to,  for  example,  external  radiation  decreases  with  decreasing  time  interval. 
However  introducing  additional  ionization,  for  example,  by  irradiation  with  a  UV  light 
source,  can  decrease  the  breakdown  voltage  of  a  pulsed  discharge.  If  the  pulse  repetition 
rate  is  high  enough  so  that  there  is  remnant  ionization  from  a  previous  discharge  pulse, 
the  pulsed  breakdown  voltage  will  also  decrease.  Therefore,  remnant  ionization  in  a 
repetitively  pulsed  DBD  or  external  ionization  will  also  influence  the  E/n  of  a  DBD. 

The  time  delay  for  the  onset  of  breakdown  after  the  application  of  the  voltage 
pulse  becomes  shorter  with  increasing  rate  of  voltage  rise  time,  and  with  the  decreasing 
delay  time,  the  voltage  at  breakdown  will  also  be  higher  for  the  faster  rise  times. 
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For  applications  using  conventional  circuit  layouts,  i.e.  when  the  energy  storage 
switch,  conducting  leads,  and  the  discharge  electrodes  are  not  all  part  of  impedance 
matched  transmission  line  configuration,  the  circuit  inductance  and  capacitance  will  have 
to  be  considered.  Very  early  investigations  of  ozonizer  discharges  have  already 
recognized  this  circuit  effect  .  The  fact  that  the  circuit  components  and  in  particular  the 
electrode  capacitance  can  have  a  considerable  influence  not  only  on  the  electrical 
characteristics  but  also  on  the  plasma  characteristics  of  the  discharge  is  the  subject  of  this 
note. 

2.  Experiment  and  electrical  measurements 

A  photograph  of  the  experimental  layout  used  in  our  measurements  is  shown  in  fig.  3.  A 
fast  high  voltage  stacked  CMOS  switch  is  used  to  apply  a  shorting  pulse  to  the  top 
electrode  of  the  two  1  x  10  cm  copper  electrodes  straddling  a  1  cm  diameter  discharge 
tube.  The  top  electrode  is  charged  to  the  high  voltage  via  a  charging  resistor.  The 
discharge  tube  is  part  of  a  vacuum  system  with  electronic  gas  flow  controllers  and  an 
automatic  pressure  controller.  The  pulse  repetition  rate  could  be  varied  from  essentially 
zero  to  about  30  kHz.  The  discharge  was  operated  at  pressures  from  30  to  200  Torr  with  a 
mixture  of  N?/Ar  (10:2)  gas.  At  the  highest  operating  pressure  range,  at  12  kV  applied 
voltage  across  the  1  cm  gap,  breakdown  required  prior  operation  at  lower  pressure  and  a 
gradual  pressure  increase  to  provide  remnant  ionization.  The  capacitance  of  the 
electrodes  was  27.5  pF,  as  detennined  from  the  electrical  characteristics. 

For  a  baseline  estimate  of  the  electrical  parameters  without  a  discharge,  the 
current  and  voltage  was  measured  at  a  pressure  of  700  Torr  and  at  the  range  of  voltages 
used  later  when  operating  at  lower  pressure  with  a  discharge.  The  measurements  were 
then  repeated  at  pressure  of  50  and  100  Torr  discharge,  as  shown  in  fig.  4.  With  a 
constant  voltage  of  12  kV  and  a  constant  pulse  repetition  frequency  of  10  kHz,  while  the 
measured  voltage  essentially  shows  no  change,  the  measured  discharge  current  adds  to 
the  displacement  current.  As  the  pressure  decreases,  the  gas  breaks  down  earlier.  This 
means  that  discharges  at  the  lower  pressures  occur  when  the  input  voltage  is  at  a  higher 
value,  so  the  discharge  operates  at  a  higher  E/n.  On  the  other  hand  the  net  current  into  the 
discharge  is  lower  for  the  lower  pressure.  When  the  operating  voltage  is  varied  while 
keeping  the  pressure  constant,  a  somewhat  different  behavior  is  observed.  Figure  5  shows 
the  discharge  currents  measured  for  different  operating  voltages  and  at  50  and  100  Torr 
pressures.  For  these  figures  the  contribution  of  the  displacement  current  for  the  respective 
operating  voltage  was  subtracted  so  that  only  the  net  current  into  the  discharge  is  shown. 
For  both  pressures,  the  current  pulse  amplitude  increases  for  increasing  operating  voltage. 
The  E/n  of  the  discharge  increases  for  increasing  operating  voltage  not  only  because  the 
operating  voltage  is  higher  but  also  because  the  onset  of  the  discharge  current  occurs 
earlier  and  at  a  higher  voltage.  (However,  the  effect  of  increasing  the  operating  voltage 
on  the  time  delay  decrease  between  the  applied  voltage  and  breakdown  becomes  smaller 
with  increasing  voltage  and  is  practically  absent  at  50  Torr  (fig. 5)).  The  difference  in 
breakdown  time  delay  and  pulse  current  amplitude  for  the  two  operating  pressures  is 
clearly  shown  in  the  inset  of  figure  5. 
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3.  Optical  measurements. 


In  order  to  investigate  the  effect  of  the  measured  electrical  characteristics  on  the  direct 
electron  impact  excitation  in  a  gas  discharge,  the  plasma  emission  intensity  of  three 
transitions  in  the  10:2  ISB/Ar  gas  mixture  were  measured.  The  three  transitions  were 
selected  for  their  difference  in  excitation  threshold  which  permits  to  estimate  the 
variation  of  the  tail  of  the  EEDF,  since  the  overlap  of  excitation  cross  sections  and  the 
EEDF  are  proportional  to  ve  fe(ve),  where  ve  is  the  velocity  of  the  electron  and  fe(ve)  is 
the  electron  energy  distribution  function.  The  net  direct  electron  impact  excitation  rates13 
are  a  product  of  both  the  electron  density  and  ve  fe(ve).  The  transitions  were:  N2  B-X 
(band  head  at  391.4  nm,  threshold  19  eV),  Ar  2pi  to  2s2  750.4  nm  (threshold  13.5  eV) 
and  N2  C-B  (threshold  1 1  eV).  The  operating  pressure  was  varied  from  30  to  200  Torr.  A 
1  m  spectrometer  with  a  Hamamatsu  R758  photomultiplier  was  used.  Figure  4  shows 
optical  intensities  (arbitrary  scale)  for  the  three  transitions  over  the  pressure  range  of  30 
to  200  Torr.  The  transition  with  the  highest  excitation  threshold,  N2+  B-X,  shows  a  steep 
drop  with  increasing  pressure  which  is  the  expected  behavior  for  a  discharge  with  a  few 
eV  mean  electron  energy  when  the  lowering  of  the  E/n  decreases  the  tail  of  the  EEDF; 
therefore  the  overlap  of  excitation  cross  sections  with  the  EEDF  also  decreases  rapidly. 
The  two  electronic  transitions  with  lower  excitation  threshold  energy,  Ar  750.4  nm  line 
and  N2  C-B,  both  show  an  initial  increase  with  increasing  pressure.  This  is  a  clear 
indication  of  the  increase  of  pulse  current  when  the  pressure  increases,  as  demonstrated  in 
fig.  4  and  5.  The  decrease  in  E/n  is  compensated,  at  least  initially,  by  the  increase  in  pulse 
current  and  power  deposition.  This  effect  is  dependent  on  the  excitation  threshold  of  the 
particular  transition.  Since  the  N2  C-B  transition  has  a  lower  threshold  than  the  Ar 
transition,  it  peaks  at  a  higher  pressure  and  it  is  less  sensitive  to  a  decreasing  E/n.  For  the 
higher  threshold  N2+  B-X  transition,  the  decreasing  discharge  E/n,  even  at  the  lowest 
pressure,  impacts  the  direct  electron  impact  excitation  rate  in  spite  of  the  increasing 
current  because  the  overlap  of  the  excitation  cross  section  and  the  tail  of  the  EEDF 
decrease  faster  than  the  electron  density.  These  emission  data  clearly  show  the  impact  of 
increasing  power  deposition,  even  when  the  E/n  decreases,  with  increasing  pressure 
between  50  and  100  Torr.  These  optical  emission  measurements  support  the  conclusions 
arrived  from  the  time  resolved  electrical  power  deposition  measurements  and  the 
estimated  E/n  values  in  a  short  pulsed  dielectric  barrier  discharge. 

4.  Conclusions 

With  fast  pulse  excitation  of  a  DBD  the  effect  of  the  displacement  current  caused  by  the 
capacitance  of  the  electrode  structure  cannot  be  neglected.  For  excitation  processes  with  a 
low  energy  threshold,  when  the  gas  pressure  is  increased  and  the  E/n  is  lowered,  the 
increase  in  net  current  into  the  discharge  can  actually  increase  the  direct  electron  impact 
excitation  rate.  The  displacement  current  limits  the  increase  of  the  discharge  current 
when  increasing  the  pulse  voltage,  since  that  also  will  increase  the  displacement  current. 
For  fast  pulse  DBD  excitation,  minimizing  the  electrode  capacitance  will  lower  the 
dielectric  component  of  the  charge  current  and  the  ohmic  losses  in  a  charging  resistor, 
therefore  increase  the  discharge  excitation  efficiency. 
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Figure  3.  Photography  of  the  DBD  setup  with  pulser. 


Upper  Electrode  Discharge 


Current 

Sensor 


High  Voltage 
Switch  w 


9 


Figure  4.  Discharge  Voltage  and  discharge  currents  for  two  pressures  for  a  N2:Ar  10:2 
gas  mixture  at  12  kV,  10  kHz.  Zero  time  has  been  arbitrarily  set  to  the  beginning  of  the 
current  pulse. (a)  50  Torr  ;  (b)  100  Torr.  Thick  line:  Total  current;  Thick  dashed  line: 
Voltage.  Thin  dashed  line:  Displacement  current  only  (700  Torr) 
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Discharge  Voltage  (kV)  Discharge  Voltage  (kV) 


Figure  5.  Net  discharge  currents  (displacement  currents  subtracted  from  measured 
current)  with  N2:  Ar  10:2  gas  mixture  at  10  kHz  and  at  different  operating  voltages,  (a): 
50  Torr,  (b):  100  Tor.  Thick  solid  line:8  kV,  Thin  dashed  line:  10  kV  (shown  vertically 
displaced).  Thin  solid  line:  12  kV  (shown  vertically  displaced)The  inset  shows  the  delay 
of  the  current  pulse  as  measured  to  the  half  amplitude  versus  the  pulse  voltage. 
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Intensity  (arbitrary  scale) 


Figure  6.  Optical  emission  intensities  for  N2  and  Ar  transitions  for  the  same  gas  mixture 
at  12  kV,  10  kHz. 
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IV.  Absolute  H  atom  density  measurement  in  a  pure  methane  pulsed 

discharge. 


In  supersonic  combustion  (scramjet)  systems  the  flow  speed  is  very  high  so  that 
mixing  and  reaction  times  are  limited  '  .  Some  flow  conditions  also  give  rise  to  very  low 
mixture  temperatures1.  For  either  condition,  additional  sources  of  ignition  or  combustion 
enhancement  are  needed.  Plasma  torch  ignitors  ’  have  been  investigated  widely,  and 
microwave  discharges3  have  been  studied,  both  as  means  of  providing  reliable  scramjet 
combustion.  However,  such  approaches  require  multi-kilowatt  electrical  power  inputs 
and  are  not  easily  volume  scaleable.  Highly  nonequilibrium  plasma  production  of 
energetic  reactive  hydrocarbon  fragments  could  have  application  in  situations  where 
more  conventional  means  of  promoting  and  sustaining  combustion  are  degraded  4’ 5. 

The  attractiveness  of  high  E/n,  where  E  is  the  electric  field  and  n  is  the  gas 
density,  pulsed  discharges  for  dissociation,  and  electronic  excitation  and  ionization  is 
clear  from  both  theoretical  ’  and  experimental  work  performed  on  corona  and 
dielectric  barrier  discharges.  Pulsed  dc  discharges  in  gaseous  hydrocarbon  mixtures  can 
provide  substantial  amounts  of  atomic  hydrogen,  as  well  as  electronically  excited 
hydrocarbon  fragment  neutrals  and  ions.  Such  discharge  products  would  be  useful  in 
promoting  combustion  in  a  wide  variety  of  fuel  to  air  mixtures. 

The  plasma  volume  and  energy  scaling  of  highly  nonequilibrium  plasma  assisted 
ignition  could  improve  overall  efficiency  of  combustion.  It  is,  therefore,  important  to 
quantify  the  radical  production  efficiency  of  a  plasma  device  which  could  be  used  to 
replace  a  conventional  spark  ignition  module.  We  have  developed  a  short  pulse  excited 
discharge  device  9  which  can  operate  at  high  reduced  electric  field  (E/n)  with  PD,  P  is  the 
gas  pressure  and  D  is  the  gap  distance,  scaling  up  to  50  Torr.cm.  Actinometric 
measurements  performed  in  this  high  reduced  electric  field  plasma  device  indicated  that 
the  CH4  dissociation  efficiency  can  be  more  than  an  order  of  magnitude  greater  than 
thennal  dissociation  process  9.  Although  those  measurements  gave  indications  that  a 
short  pulse  discharge  can  produce  radicals  with  higher  efficiency,  compared  to  the 
thennal  dissociation,  those  measurements  were  performed  at  low  current  density  (2 
Amp/cnr)  and  low  energy  density  (0.7mJ/cm2).  In  this  paper  we  report  absolute  H  atom 
ground  state  density  measurements  in  a  short  pulse  methane  discharge  by  two-photon 
allowed  LIF  (TALIF)  technique  at  current  densities  exceeding  20  Amp/cnr.  The  H  atom 
TALIF  intensity  has  been  calibrated  by  comparing  TALIF  intensity  from  Kr  atom  ground 
state  10  to  obtain  absolute  H  atom  density. 

The  goals  of  this  study  are:  1)  construct  a  high  E/n  discharge  system  capable  of 
producing  discharges  over  a  substantial  range  of  input  energy  density,  and  2)  conduct 
spectroscopic  measurements  of  the  absolute  H-atom  production  to  estimate  plasma 
dissociation  efficiency  of  CH4  in  high  E/n  discharges.  Pulsed  gas  discharges  in  CH4 
produce  all  possible  methane  fragments  (CH3,  CH2,  CH,  C  and  H)  and  their  ions  in 
varying  ratios  depending  upon  the  discharge  E/n.  In  the  present  study,  we  provide  an 
accurate  estimate  of  the  direct  electron  impact  dissociation  efficiency  of  CH4  in  short 
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pulse  excited,  high  E/n  discharge,  with  nearly  two  orders  of  magnitude  higher  energy 
density  compared  to  the  previous  measurement 9 ,  using  absolute  H  atom  number  density 
measurement  by  calibration  of  H  atom  TALIF  signal  with  near  isoenergetic  Kr  ground 
state  atom  TALIF  measurement 10. 

A  solid-state,  high-voltage,  high-speed  switch  (15  kV,  30  A)  provided  a  fast  rise 
time  voltage  pulse  and  was  controlled  by  timing  control  electronics.  As  shown  in  Figure 
7,  the  voltage  pulse  was  delivered  (rise  time  <15  ns)  via  a  microstrip  transmission  line  to 
an  end-fed  coaxial  section. 

Four  fused  silica  windows  provided  optical  access.  The  discharge  gap  was  viewed 
by  a  gated  photomultiplier  tube  (PMT)  fitted  with  an  optical  relay  lens  TL  and  preceded 
by  F  a  long  pass  and  a  neutral  density  filter.  The  gated  PMT  allows  photon  detection  after 
the  plasma  voltage  pulse  was  turned  off  and  the  plasma  emission  had  decayed.  The  laser 
was  fired  185  ns  after  the  voltage  pulse  was  turned  off.  The  pulse  discharge  and  the 
TALIF  data  acquisition  timing  was  controlled  by  a  delay  generator  pulser  and  t=0  on  all 
the  data  plots  correspond  to  the  start  of  the  voltage  pulse. 

The  615  nm  pulsed  dye  laser  output  tripled  by  a  combination  of  KDP  and  BBO 
crystals  and  the  TALIF  emission  was  measured  with  a  gated  PMT.  The  optical  path, 
including  attenuation  and  collection  solid  angle,  was  kept  fixed  for  all  measurements.  The 
laser  beam  was  focused  near  the  middle  of  the  discharge  volume  and  it  was  located  2  mm 
below  the  cathode  surface  for  optimum  TALIF  signal.  Absolute  calibration  used 
successive  measurements  of  TALIF  signals  at  two  wavelengths;  205.08  nm  probed  the 
ground  state  of  hydrogen  atom  produced  by  the  pulsed  discharge  in  methane  and  the 
204.13  nm  TALIF  measured  the  calibration  signal  from  the  ground  state  of  krypton 
atom.10  These  two  probe  wavelengths  are  close  enough  to  ensure  the  input  optics  spectral 
throughput  parameters  remain  constant.  The  205.08  nm  wavelength  two-photon  excites 
population  of  hydrogen  atoms  from  the  ground  state,  Is,  to  3d  and  3s  levels.  The  TALIF 
signal  is  detected  by  photon  emission  at  656.3  nm  (Ha).  The  other  probe  wavelength 
excites  a  krypton  transition  from  the  ground  state,  4p6,  to  5p  levels,  using  two  photons  at 
204.13  nm  and  the  emission  from  5p  to  5s7  was  detected  at  826.3  nm.  For  ease  of 
calibration  the  two  TALIF  signals  were  set  approximately  equal  amplitude  by  adjusting 
the  pressure  of  Kr.  A  pyroelectric  detector  measurement  was  conducted  with  the  series  of 
TALIF  measurements  to  assure  that  probe  laser  energy  had  not  changed  more  than  5%, 
which  limits  the  TALIF  measurement  error  to  10%  due  to  the  laser  intensity  variation. 

All  measurements  have  been  perfonned  in  pure  CH4  discharge  at  20  Torr  pressure 
with  5  mm  gap  between  the  1  cm  diameter  electrodes.  The  applied  voltage  was  varied 
from  4  kV  up  to  6  kV  with  480  ns  pulse  duration.  Gas  flow  and  pressure  was  controlled 
by  a  mass  flow  controller  loop  using  an  exhaust  valve  controller  and  a  mass  flow  readout, 
and  a  Baratron  pressure  transducer.  Throughout  this  study  the  gas  flow  rate  was  fixed  at 
30  standard  cubic  centimeters  per  minute.  The  CH4  gas  was  ultra  high  purity  grade 
(99.999%). 
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The  pulsed  discharge  voltage,  and  current  deposited  in  the  discharge  are  shown  in 
Figures  8a,  and  8b,  respectively.  The  discharge  voltage  measurement  of  high  speed 
pulses  (<15  ns)  show  the  gap  voltage  initially  spikes,  due  to  pre  breakdown  over  volting, 
to  a  high  value.  At  6  kV  applied  voltage,  the  20  Torr  pure  methane  discharge  runs  at 
current  density,  j,  up  to  25  Amp/cm  .  The  j/P  scaling  of  the  cathode  sheath  voltage 
drop1 1  for  pure  methane  discharge  is  similar  to  that  of  FL  discharge;  the  cathode  fall 
voltage  can  be  700-750  V  for  these  conditions.  (Measurement  of  the  discharge  voltage 
using  a  variable  gap  discharge  confirms  the  estimated  sheath  voltage  drop).  About  20  to 
50  ns  after  breakdown,  depending  on  drive  voltage,  a  lower  plateau  gap  voltage  and  the 
discharge  impedance  is  established  for  the  measured  480  ns  voltage  pulse  duration.  The 
discharge  operating  voltage  is  dependent  on  the  initial  applied  voltage;  it  varies  from 
1 100  V  at  6  kV  (sheath  voltage  drop  700  V)  down  to  900  V  at  4  kV  (sheath  voltage  drop 
490  V).  The  corresponding  volume  average  discharge  E/n,  assuming  gas  temperature  is 
400  K,  can  vary  from  85  Td  to  90  Td  (lTd  =  10~17  Volt  cm'2).  The  peak  E/n,  near 
breakdown,  can  exceed  300  Td. 


Pressure  dependent  krypton  TALIF  calibration  signal  were  measured  and  was 
found  to  be  linear  with  pressure.  Calibration  transfer  was  made  via  a  measurement  at 
0.3 18  Torr,  where  Kr  signal  amplitude  was  comparable  to  the  H  atom  TALIF.  A 
representative  set  of  this  series  is  shown  in  Figure  9  for  three  pressures  0.098,  0.191,  and 
0.318  Torr.  Immediately  following  the  Kr  atom  TALIF  measurement,  ground  state  H 
atom  production  in  CFfi  discharge  was  measured. 


The  H  atom  TALIF  intensity  has  been  calibrated  by  comparing  TALIF  intensity 
from  Kr  atom  ground  state  10  to  obtain  absolute  H  atom  density  using  the  equation  given 
below. 


Hh 


„  a2i(.Kr)  cr<2)(Kr)  SPMT(H)f  I £(ifr)Y 
S  PMT  (Kr){lL(H)J 


(1) 


where  nH  and  nKr  are  the  absolute  number  density  of  the  H-  and  Kr-atoms  respectively,  k 
is  a  transmission  and  sensitivity  factor  for  the  optical  set-up  and  detector  at  the  two 
wavelengths,  a23  is  the  optical  branching  ratio  (a23=  A23  /  X  A2,,  A23  is  the  spontaneous 
emission  coefficient  and  XA2]  is  the  sum  of  spontaneous  emission  and  collisional 
quenching  coefficients),  a<2)  is  the  two-photon  cross  section,  Spmt  is  the  measured  TALIF 
signal,  and  II  is  the  intensity  of  the  laser  radiation.  The  ratio  of  the  two-photon  cross- 
sections,  Kr  atom  optical  branching  ratio  and  collisional  quenching  rate  of  H  atom  3d  and 
3s  states  by  CH4  are  given  in  reference  10.  Under  our  measurement  conditions,  the  value 
of  k  depends  on  the  ratio  of  the  PMT  response  sensitivity  at  656.3  nm  and  826.3  nm, 
which  is  obtained  from  the  manufacturer’s  data  sheet.  Thus  from  a  measurement  of  both 
the  H-atom  and  Kr-atom  TALIF  signals,  the  absolute  number  density  of  H-atoms  can  be 
calibrated  with  a  known  number  density  of  ground  state  Kr  atom.  Figure  10  shows  H 
atom  TALIF  signals,  calibrated  in  absolute  H  atom  number  density  by  comparing  with  Kr 
atom  TALIF  signal,  obtained  with  three  different  applied  drive  voltages  from  4  kV  to  6 
kV,  with  480  ns  pulse  duration.  The  total  input  energy  for  these  discharges  were  6.5,  8.7, 

1 1.4  mJ,  respectively.  For  the  corresponding  discharge  conditions,  the  measured  H  atom 
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density  ranged  from  ~  9.8x  1015  per  cm3  at  4  kV  to  2.2*  1016  per  cm3  at  6  kV  (TALIF  pulse 
peak-to-area  ratio  is  essentially  constant  so  the  peak  values  are  used),  i.e.,  nearly  4.5%  of 
CH4  was  dissociated  at  6  kV  applied  voltage,  assuming  400  K  gas  temperature  at  20  Torr, 
across  the  5  mm  gap  discharge.  Based  on  the  energy  deposited  in  the  discharge  volume 
near  the  cathode,  an  estimate  of  the  direct  electron  impact  dissociation  efficiency  of  CH4 
+e  — >  CH3  +  H  +e  can  be  obtained  from  the  measured  H  atom  yield  in  this  pulsed 
discharge.  The  product  of  the  direct  electron  impact  dissociation  rate  coefficient  and 
electron  density  in  plasma  is  proportional  to  the  fraction  of  discharge  energy  deposited  to 
that  process  multiplied  by  the  energy  deposited  in  the  gas,  i.e.: 
kTn™  *Y(E. In)  PCH<  (2) 

where  Y(E  /  n )  is  the  fractional  energy  deposited  into  electron  impact  production  of 
CH3+H(n=l)  as  a  function  of  E/n  in  methane  discharges,  k^H4  is  the  electron  impact 
dissociation  rate,  n^HA  is  the  electron  density,  and  PCH4  is  total  discharge  energy.  For  this 

discharge  operating  condition  with  700  V  sheath  voltage  drop  along  with  about  80  V/cm 
plasma  electric  field  in  the  5  mm  gap  discharge,  70  to  75%  of  the  energy  is  deposited  in 
the  discharge  volume  within  2  mm  from  the  cathode  surface.  For  the  6  kV  applied 
voltage,  discharge  with  1 1 .4  mJ  total  discharge  energy,  this  would  correspond  to  the 
discharge  energy  density  ~  70  mJ/cm3.  Assuming  an  average  CH4  dissociation  energy 
cost  of  10  eV,  the  total  energy  required  to  produce  H  atom  density  of  2.2x  10  per  cm 
would  be  35  mJ/cm  (direct  electron  impact  dissociation  threshold  “  of  CH4+  e  — >  CH3  + 
H  +  e  is  8.5  eV,  and  at  the  threshold  energy  the  dissociation  rate  is  very  small.  The  higher 
dissociation  rate,  k^H4 ,  with  average  electron  energy  cost  of  10  eV  would  be  more 
appropriate  for  this  short  pulse  discharge).  This  would  imply  about  50%  efficiency  of  the 
direct  electron  dissociation  of  CH4  in  this  6  kV  short  pulse  discharge.  Similarly  for  the  4 
and  5  kV  discharges,  the  corresponding  energy  densities  are  ~  40  and  53  mJ/cm3, 
respectively.  The  energy  required  to  produce  9.8x10  ~  per  cm  H  atom  density,  at  4  kV, 
would  be  15.7mJ/cm3  and  the  CH4  dissociation  efficiency  is  ~  40%.  For  the  5  kV 
discharge,  the  energy  required  to  produce  H  atom  density  of  1 .44x10  per  cm  would  be 
23  mJ/cm3  and  the  dissociation  efficiency  is  nearly  the  same  as  the  4  kV  discharge.  The 
measurable  difference  in  the  CH4  dissociation  efficiency  at  higher  applied  voltage  may  be 
caused  by  the  increase  in  the  discharge  energy  deposition  in  the  sheath,  where  nonlocal 
electron  impact  excitation  can  dissociate  CH4  at  higher  efficiency  than  in  the  volume 
plasma.  Note  that  volume  plasma  E/n  is  essentially  constant  for  all  the  applied  voltages. 
Although  the  importance  of  initial  energy  deposition  near  the  breakdown  under  higher 
E/n  can  not  be  directly  quantified  from  these  measurements,  these  measurements  clearly 
show  the  advantage  of  operation  at  higher  applied  voltage  even  for  a  relatively  long  pulse 
duration  discharge  compared  to  the  initial  high  peak  power  pulse  duration.  In  summary, 
we  have  used  absolute  H  atom  density  measurements  by  TALIF  to  demonstrate  high 
efficiency  operation  of  short  pulse  excited  discharge  for  CH4  dissociation,  H  atom  and 
CH3  radical  production  which  can  be  used  for  ignition  and  combustion  enhancement  for 
high  speed  propulsion. 
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FIGURES 


Figure  7.  An  experimental  schematic;  a)  laser  beam  path  through  a  focusing  lens  L  and 
quartz  windows  W,  b)  a  housed  PMT  with  a  view  of  the  gap  orthogonal  to  the  beam 
through  transfer  lenses  TL  and  long-pass/neutral  density  filters  F  and  c)  a  microstrip 
transmission  line  MT  feeding  the  gap  on  center  of  a  symmetrical  pair  of  cylindrical 
section  copper  conductors,  the  conventional  current  I  is  shown  with  arrows. 
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Figure  8a.  Discharge  gap  voltage  for  480  ns  pulse  width  4  and  6  kV  applied  voltages 
across  a  5  mm  gap,  1  cm  diameter  electrodes,  at  20  Torr  CH4. 


Figure  8b.  Discharge  current  pulses  for  480  ns  pulse  width  4,  5,  and  6  kV  applied 
voltages,  across  a  5  mm  gap,  1  cm  diameter  electrodes,  at  20  Torr  CH4. 
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Figure  10.  Krypton  atom  calibration  TALIF  measurements  for  Kr  gas  pressures  0.098, 
0.3191,  and  0.318  Torr. 


Figure  11.  TALIF  measurement  of  absolute  H  atom  density  in  20  Torr,  480  ns  pulse 
width,  CH4  discharge  with  4,  5,  and  6  kV  applied  voltages.  Absolute  H  atom  densities  are 
obtained  from  the  Kr  atom  TALIF  calibration,  yielding  H  atom  density  of  9.8 x  1 015  per 
cm3  at  4  kV,  1.44><1016per  cm3  at  5  kV,  and  2.2><1016per  cm3  at  6  kV. 
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V.  Evidence  of  enhanced  electronic  excitation  at  the  shock  front  of  low 
Mach  number  shocks  in  non-equilibrium  plasmas 


The  propagation  and  dispersion  of  weak  shock  waves  in  non-equilibrium  plasmas 
have  recently  received  considerable  attention  '"6.  Optical  and  pressure  probe 
measurements  show  spreading  of  the  shock  front  and  an  increase  in  shock  velocity  in  the 
plasma.  Some  investigators  relate  these  effects  to  axial  or  radial  thermal  gradients4’5  in 
the  plasma,  others  invoke  effects  caused  by  a  shock  wave  induced  electrostatic 
interactions  ’  in  the  plasma.  Depending  on  the  boundary  conditions  of  the  experiments, 
these  interpretations  may  be  applicable.  We  have  recently  investigated  6  the  interaction  of 
the  shock  wave  with  the  N2  plasma  in  a  DC  glow  discharge  and  have  shown  that  the 
propagating  shock  wave  causes  large  local  changes  in  the  electric  field  in  the  immediate 
vicinity  of  the  shock  front  as  well  as  global  changes  in  discharge  voltage  and  current.  The 
magnitudes  of  these  effects  were  shown  to  be  dependent  on  discharge  polarity.  Spatially 
and  temporally  resolved  measurements  of  the  optical  emission  from  the  C  1TU-  B  ng 
transition  in  nitrogen  showed  sharp  decreases  when  the  shock  wave  passed  through  the 
observation  point.  Fig.  12  shows  a  series  of  ICCD  camera  images  taken  with  a  1.2  ps 
exposure  time  and  a  blue  filter,  which  eliminates  the  emission  of  the  nitrogen 
recombination  afterglow.  The  ICCD  pictures  were  made  on  a  5  cm  diameter,  5  Torr,  60 
mA  discharge. 

As  the  shock  wave  sweeps  through  the  positive  column,  the  discharge  becomes 
dark  behind  the  shock.  This  result  indicates  a  significant  decrease  in  the  reduced  electric 
field  E/n,  where  E  is  the  electric  field  and  n  is  the  neutral  density,  behind  the  shock  even 
though  there  is  continuity  of  current.  Note  also  that  the  emission  intensities  near  both  the 
cathode  and  the  anode  show  no  discernible  change  in  intensity,  confirming  that  under 
these  operating  parameters  the  discharge  current  changes  little.  These  observations  are 
correlated  with  an  observed  sharp  decrease  of  the  reduced  electric  field  at  the  shock 
arrival  time  as  measured  by  electrical  probes  6.  The  decrease  in  optical  emission  to 
practically  zero  intensity  behind  the  shock  front  indicates  a  reduction  of  electron  impact 
excitation  and  therefore  also  of  ionization.  As  already  observed  from  the  ICCD  images  of 
the  cathode  and  anode  glow,  the  discharge  current  decreased  only  moderately  during  the 
time  when  the  shock  passed  through  the  length  of  the  positive  column  discharge  tube  for 
anode  to  cathode  propagation.  It  even  showed  a  small  temporary  increase  for  shock 
propagation  in  the  direction  from  cathode  to  anode.  The  decrease  of  electronic  excitation 
and  ionization  could  be  explained  by  the  large  increase  of  neutral  density  behind  the 
shock  front  and  therefore  a  decrease  of  E/n.  However,  the  decrease  of  E/n  with  the 
increase  in  neutral  density  will  also  reduce  the  electron  drift  velocity.  Therefore, 
discharge  current  continuity  requires  an  increase  of  electron  density,  compensating  drift 
velocity  decrease,  by  a  factor  of  2  for  the  Mach  1.7  shock  waves.  The  observed  sharp 
drop  in  the  local  electric  field  propagating  with  the  shock  front  indicated  a  local  region  of 
enhanced  conductivity  6.  This  abrupt  change  in  electric  field  (requires  net  space  charge  to 
satisfy  Poisson’s  equation)  suggested  the  formation  of  a  triple-  or  quadruple  space  charge 
layer  connected  with  the  shock  front.  The  high  electric  fields  in  this  very  localized  space 
charge  region,  ahead  of  the  shock  front,  can  then  lead  to  local  (dimension  determined  by 
the  plasma  Debye  length)  excitation  and  ionization  that  can  provide  the  excess  electron 
density.  The  calculated  electron  density  decay  in  an  equivalent  unperturbed  positive 
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column  discharge  was  estimated  to  be  at  most  20%  during  a  time  equal  to  the  shock 
traversal  time. 

The  optical  and  electrical  measurements  gave  a  self-consistent  picture  of  the 
plasma-shock  wave  interactions  and  the  new  measurements  reported  here  provide 
evidence  for  the  presence  of  strong  excitation  (and  ionization)  enhancement  in  a  very 
narrow  region  propagating  with  the  shock  wave. 

Under  the  conditions  of  the  discharge  (5  Torr,  50- 100mA),  the  upstream  neutral  gas  mean 
free  path  (temperature  corrected)  is  0.005cm  and  the  estimated  Debye  length  for  the 
largest  current  is  0.01cm  6.  Assuming  that  the  increase  in  optical  emission  will  occur  over 
the  Debye  length  connected  with  the  plasma  conditions  of  the  passing  shock  front,  it  will 
have  a  width  of  less  than  0.3  mm.  The  spatial  resolution  requires  that  the  full  width  of  the 
discharge  be  within  the  depth  of  field  of  the  photomultiplier  optics.  Taking  a  discharge 
diameter  of  4  cm,  on  the  outer  limits  of  the  diameter  of  the  discharge,  the  optics  used  (f 
number  30,  90  mm  focal  length)  will  allow  a  resolution  of  only  1.7  mm  (or  blur  disc 
diameter)7.  A  resolution  of  0. 1  mm  is  possible  over  a  depth  of  field  of  2.4  mm  only.  The 
combined  effects  of  short  radiation  enhancement  time  and  insufficient  spatial  resolution 
will  broaden  the  pulse  width  of  a  shock  induced  enhanced  optical  radiation  pulse  and 
reduce  its  amplitude,  making  its  detection  difficult.  These  limitations  may  have  been 
responsible  for  our  inability  to  detect  the  expected  electron  impact  excitation 
enhancement  in  C  to  B  emission. 

In  order  to  improve  our  plasma  emission  detection  sensitivity,  we  have  measured 
the  2-0  vibrational  transition  of  the  B  ng  -  A  Zu  first  positive  band,  which  has  a 
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radiative  lifetime  of  8.46  ps  .  This  longer  lifetime  along  with  the  lower  excitation  energy 
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level  compared  to  the  C  nu  state  was  expected  to  give  a  better  chance  to  observe  the 
enhanced  emission  peak  at  the  shock  front.  The  B  state  can  also  be  populated  by 
cascading  from  three  other  states.  The  experimental  setup,  except  for  the  spectral  filter 
and  a  3  cm  diameter  discharge  tube  with  electrodes  conformal  with  the  glass  wall,  was 
identical  to  the  one  used  previously  6 .  This  discharge  geometry  permits  shock  wave 
propagation  measurements  with  the  discharge  polarity  reversed  with  minimal 
modification  of  the  test  setup.  The  tube  diameter  was  3  cm  diameter  vs.  5  cm  in  the 
previously  reported  measurements  6,  and  at  3  Torr,  the  Mach  number  is  about  the  same 
than  in  the  previous  measurements  (Mach  1.76)  using  100  Joules  spark  gap  energy  input. 
The  narrower  tube  diameter  and  lower  pressure  also  results  in  increased  electron  losses  to 
the  wall  in  this  diffusion-dominated  discharge  so  the  electron  mean  energy  is  slightly 
higher  than  in  the  discharge  of  our  previous  report6.  The  optical  observation  point  was 
118  mm  from  the  electrode  closer  to  the  spark  gap  and  the  distance  between  the 
electrodes  was  180  mm.  Fig.  13a  shows  the  discharge  voltage,  current  and  775.4  mn 
emission  for  a  shock  wave  propagation  from  cathode  to  anode,  fig.  13b  from  anode  to 
cathode,  in  both  cases  for  a  DC  current  of  20  mA  at  3  Torr.  The  laser  photo  deflection 
markings  of  the  shock  front  arrivals6  in  the  775.4  mn  emission  diagrams  are  for  the 
physical  location  of  the  photomultiplier.  The  time  period  during  which  the  shock  wave 
traverses  the  discharge  is  indicated  by  the  vertical  dashed  lines.  As  observed  previously, 
the  discharge  voltage  decreases  and  the  current  increases  by  10%  for  about  100  ps  for 
cathode  to  anode  shock  propagation  direction,  for  the  opposite  propagation  direction  the 
voltage  increases  and  the  current  decreases  immediately.  Again  as  previously  observed, 
locally  at  the  photomultiplier  observation  point  the  light  intensity  increases  when  the 
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shock  wave  enters  the  discharge  for  cathode  to  anode  propagation  and  decreases  for  the 
opposite  propagation  direction.  The  electrical  as  well  as  the  optical  effects,  before  arrival 
of  the  shock  at  the  observation  point,  are  caused  by  the  global  electrical  circuit  response 
to  the  propagating  shock  wave,  as  described  previously  6.  The  global  voltage  and  current 
perturbations  in  this  smaller  diameter,  lower  pressure  discharge  are  larger  than  for  the 
previous  measurements  6.  A  sharp  and  pronounced  peak  of  775.4  mn  emission  is 
observed  when  the  shock  wave  passes  through  the  observation  region  for  both 
propagation  directions.  The  emission  peaks  have  a  rise  time  of  2  ps  and  a  decay  time  of 
about  40  ps.  This  long  decay  is  in  contrast  to  the  337.1  mn  emission  observed  previously, 
which  decays  within  less  than  2  ps  to  practically  zero  intensity  in  conjunction  with  the 
fast  increase  in  neutral  density  behind  the  shock.  The  long  decay  time  is  probably  caused 
by  the  heavy  particle  collision  induced  excitation9  of  the  B  state  due  to  the  reaction 
N2( A3XL1+)+N2(X 1  Eg+,v>5)  ->  N2(B3ng)+N2(X1Eg+).  (Both  the  A3IL1+  state  and  the 
vibrationally  excited  ground  state  have  long  life  times). 

The  rapid  rise  in  intensities  of  the  775.4  mn  emission  is  an  indicator  of  a  large 
increase  in  electronic  excitation  rate  connected  with  the  shock  wave  induced  very  local 
increase  in  E/n  value  (or  electronic  energy  distribution  function).  Due  to  the  limitations  of 
the  optical  system  mentioned  above,  the  temporal  width  of  the  rising  portion  of  the 
measured  775.4  mn  radiation  pulses,  which  corresponds  to  the  apparent  duration  of  the 
shock  wave  generated  local  increase  in  E/n  pulses,  is  much  longer  than  the  actual  width 
of  these  radiation  pulses  and  also  the  duration  of  the  increased  electronic  excitation.  For 
the  same  reason,  the  amplitudes  of  the  measured  pulses  are  much  smaller  than  the  actual 
radiation  amplitudes.  For  this  positive  column  discharge  6,  which  runs  at  E/n  «  45  Td  (1 
Td=  lxlO'17V  cm'1),  the  observed  B3ng-A3EU4  intensity  jumps  by  at  least  20% 
compared  to  the  steady  state  in  2  ps.  This  requires  a  rate  of  change  of  the  electron  impact 
excitation  by  a  factor  of  105  indicating  that  the  E/n  increases  10  from  45  Td  to  >  240Td. 

The  decay  of  the  775.4  mn  radiation  is  strictly  due  to  the  population  transfer  from  long- 
lived  states  after  the  short  period  of  direct  electron  impact  excitation  enhancement;  no 
such  long  decay  was  observed  with  the  337. 1  mn  radiation  since  the  C  state  is  populated 
primarily  by  direct  electron  impact. 

In  fig.  14  the  emission  signals  from  the  C  to  B  and  B  to  A  states  are  compared. 

The  photomultiplier  output  traces  are  aligned  such  that  the  signal  amplitudes  before 
arrival  of  the  shock  wave  are  about  equal  on  the  plot.  The  rise  of  the  775.4  mn  signal 
corresponds  to  the  arrival  of  the  shock  wave,  and  the  decay  of  the  337. 1  mn  signal  to  the 
departure  of  the  shock  wave  from  the  observation  point.  Due  to  the  limitations  in  spatial 
resolution  discussed  above,  the  half  width  between  the  two  signals  is  about  2  ps,  which 
would  correspond  to  a  width  of  the  shock  wave  of «  1  mm  at  a  velocity  of  Mach  1.76.  As 
discussed  above,  the  actual  width  is  probably  more  than  an  order  of  magnitude  smaller. 

The  decay  of  the  discharge  current  follows  an  exponential  function.  While  this  dependence 
usually  indicates  ambipolar  diffusion  decay  to  the  wall,  the  decay  time  constant  is  smaller 
than  expected  from  a  plasma  with  an  estimated  electron  temperature  of  about  0.3  eV  at  the 
elevated  pressure  after  the  shock.  We  may  assume  that  some  of  the  charge  carriers  leave 
through  the  ends  of  the  discharge  since  in  this  experiment,  unlike  in  a  conventional  pulsed 
afterglow,  there  is  still  an  applied  electric  field.  However  the  E/n  is  much  reduced  and  too 
low  for  efficient  ionization.  Therefore,  even  with  reduced  drift  velocity,  the  current  to  the 
discharge  electrodes  still  will  cause  additional  plasma  loss. 
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The  observed  B  to  A  emission  enhancements  show  that  the  postulated  6,  triple  or 
quadruple  space-charge  layer  connected  with  a  shock  wave  propagating  in  non-equilibrium 
plasmas  can  induce  strong  local  electronic  excitation.  We  are  presently  investigating  the 
consequence  of  the  increase  in  local  E/n  on  ionization  enhancement,  and  also  the  local  Joule 
heating  that  could  be  generated,  at  or  near  the  shock  front,  by  such  an  ionization 
enhancement  produced  by  the  shock  induced  space-charge  layer  fonnation  in  a 
nonequilibrium  plasma. 
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Figure  12:  ICCD  camera  pictures  of  the  positive  column  light  emission  response  of  a  N2 
discharge  (5  Torr,  60  mA  discharge  current,  1.2ps  exposure  time,  blue  fdter,  tube 
diameter  5  cm)  to  a  shock  wave  entering  from  the  lhs  of  the  photographs.  The 
progression  of  the  shock  wave  from  cathode  to  anode  is  approximately  indicated  by  the 
dashed  line. 
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Figure  13a:  Discharge  voltage,  current  and  775.4  mn  emission  in  a  3  Torr,  20  inA 
discharge  in  N2  for  a  shock  propagating  from  cathode  to  anode.  A  laser  photodeflection 
signal  for  the  location  of  the  photomultiplier  is  shown  in  the  775.4  mn  emission  diagram. 
The  time  during  which  the  shock  traverses  the  discharge  is  indicated.  (Photomultiplier 
intensity  increases  downwards). 
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Figure  13b:  Same  as  figure  2a  but  for  anode  to  cathode  shock  propagation. 
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Figure  14:  Comparison  of  337.1  nm  and  775.4  nm  emission  signals  for  the  same 
conditions  as  in  figure  2a.  The  scales  for  each  signal  are  adjusted  such  that  the  signal 
amplitudes  before  arrival  of  the  shock  are  equal.  The  775.4  emission  increase  correlates 
with  the  arrival  of  the  shock  wave  at  the  observation  point,  the  337.1  emission  decrease 
with  the  departure.  The  vertical  markers  indicate  the  time  difference  between  the 
enhanced  excitation  (B-A),  increased  E/n  signal  and  the  signal  decrease  (C-B)  caused  by 
the  density  jump. 
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VI.  Methane  Dissociation  in  Pulsed  DC  Discharges  at  High  E/N 


Nomenclature 


E  -  electric  field  magnitude 
N  -  gas  number  density 

E/N  -  reduced  electric  field  measured  in  Townsends  (Td) 

(note:  1  Td  =  10'17  V  cm2  =  0.323  V/Torr  cm  at  300  K) 

XY(state  label  multiplicity  parity)  -  generic  molecular  electronic  state  designation 
v  -  change  in  vibrational  quantum  number 
p  -  pressure 
d  -  discharge  gap 
n  -  principal  quantum  number 
gas  -  hydrogen  or  methane 

X(E/N)  -  fraction  of  discharge  power  deposited  into  dissociation  as  a  function  of  E/N 
Pgas  -  discharge  power  deposited  into  a  given  gas 
De  -  energy  required  for  dissociation 
V  -  discharge  volume 
t  -  time 

I  Ha  '  intensity  of  H  emission  from  a  given  discharge  gas 
C  -  constant  containing  all  detection  system  optical  and  geometrical  parameters 
h  .  -  photon  energy  for  H(n=3  to  2)  or  H  emission 

An,  Asj  -  Einstein  coefficient  for  spontaneous  emission  for  the  indicated  transition 
kfs  -  electron  impact  dissociative  excitation  rate  coefficient 

nfs  -  electron  density  in  a  given  discharge  gas 

kgS  -  quenching  coefficient  for  H(n=3)  atoms  in  a  given  gas 

[Hi]  -  concentration  of  molecular  hydrogen 
[CH4]  -  concentration  of  methane 

Ef™  (DE)  -  threshold  for  dissociative  excitation  to  H(n=3)  for  a  given  discharge  gas 
Ef™  ( D )  -  threshold  for  dissociation  of  a  given  discharge  gas 

k[‘‘s  -  rate  coefficient  for  electron  impact  excitation  to  H(n=3)  for  a  given  discharge  gas 
[H] gaS  -  concentration  of  H  atoms  in  a  given  discharge  gas 

Y(E/N),  Y'(E/N)  -  fraction  of  discharge  power  deposited  into  electron  impact  excitation  as 
a  function  of  E/N 


Introduction 

Atomic  hydrogen,  hydrocarbon  fragments  and  their  radicals  are  important  species 
for  initiating  ignition  and  sustaining  stable  flame  propagation.  Molecular  gas  discharges 
generally  give  rise  to  a  complex  set  of  interactions  between  atoms,  internally  excited 
molecules  and  radicals.  Globally,  discharge  power  is  deposited  via  inelastic  electron 
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collisions  into  molecular  species  as  rotational,  vibrational  or  electronic  excitation,  and  via 
dissociation  or  ionization.  Clearly  the  various  channels  compete  for  discharge  energy, 
with  the  fractional  power  transferred  to  each  being  a  strong  function  of  E/N  (where  E  is 
the  electric  field  and  N  is  the  gas  density).  At  low  E/N  rotational  and  vibrational 
excitation  typically  dominates.  As  the  discharge  E/N  is  increased,  dissociation,  electronic 
excitation  and  ionization  collectively  share  the  bulk  of  the  input  discharge  power  ’  . 
Highly  nonequilibrium  plasma  production  of  energetic  reactive  hydrocarbon  fragments 
could  have  application  in  situations  where  more  conventional  means  of  promoting  and 
sustaining  combustion  are  degraded. 

One  such  situation  is  found  in  supersonic  combustion  ramjet  (Scramjet) 
propulsion  systems  where  the  flow  speed  is  so  high  that  mixing  and  reaction  times  are 
limited  '  .  Some  flow  conditions  also  give  rise  to  very  low  mixture  temperatures  .  For 
either  condition,  additional  sources  of  ignition  or  combustion  enhancement  are  needed. 
Plasma  torch  ignitors  have  been  investigated  widely,  and  microwave  discharges  have 
been  studied,  both  as  means  of  providing  reliable  Scramjet  combustion.  However,  such 
approaches  require  multi-kilowatt  electrical  power  inputs  and  are  not  easily  volume 
scaleable. 

The  efficacy  of  the  plasma  torch  approach  is  attributed  to  the  introduction  of  atomic 
hydrogen  into  the  flow.  A  plasma  torch,  typically  using  a  flowing  Ar-H2  mixture, 
operates  at  low  E/N  and  acts  to  thermally  dissociate  hydrogen  in  a  low  E/N  arc  discharge. 
Hence  the  energy  requirement  is  high.  Pulsed  dc  glow  discharges  with  high  E/N  in  pure 
H2  and  H2-N2  mixtures1,2  have  demonstrated  efficient  dissociation  of  hydrogen  with 
substantially  lower  power  requirements.  The  attractiveness  of  high  E/N  pulsed 
discharges  for  fostering  dissociation,  electronic  excitation  and  ionization  is  clear.  Pulsed 
dc  discharges  in  gaseous  hydrocarbon  mixtures  similarly  provide  substantial  atomic 
hydrogen,  as  well  as  electronically  excited  hydrocarbon  fragment  neutrals  and  ions.  Such 
discharge  products  would  be  useful  in  promoting  combustion  in  a  wide  variety  of  fuel  to 
air  mixtures.  Additionally,  control  of  discharge  parameters  offers  the  possibility  of 
varying  the  yields  of  the  dominant  reaction  products  in  response  to  changing  conditions 
in  an  application. 

Discharges  in  methane  and  methane-bearing  mixtures  have  been  the  focus  of  a 
number  of  recent  studies  with  a  variety  of  potential  applications.  For  example,  pulsed 
discharges  in  methane  have  been  shown  to  provide  efficient  conversion  to  acetylene  and 
hydrogen.  ’  Similar  discharges  have  been  employed  for  plasma  carburizing  of  surfaces 
and  for  the  production  carbon-based  nanoparticles10,  n. 

The  primary  goals  of  this  study  were:  1)  construct  a  discharge  system  capable  of 
producing  unifonn  methane  discharges  over  a  substantial  range  of  high  E/N  values,  and 
2)  conduct  imaging  and  spectroscopic  measurements  of  the  resulting  plasma  emissions 
sufficient  to  document  plasma  uniformity,  and  provide  an  estimate  of  H-atom  production. 
Pulsed  gas  discharges  in  CH4  produce  all  possible  methane  fragments  (CH3,  CH2,  CH,  C 
and  H)  and  their  ions  in  varying  ratios  depending  upon  the  discharge  E/N.  Additionally, 
some  carbon  clustering  occurs9'11.  Of  these,  excited  neutrals  of  atomic  hydrogen,  CH  and 
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C2  were  chosen  for  spectroscopic  observation  due  to  their  easily  detected  visible 
emissions  (H  at  656  nm,  the  CH  (A2  -X2  )  band  at  43 1  nm,  and  the  C2(A3  g-X3  u. 

v  =  0)  Swan  bands  near  516  nm  ’  ).  Comparison  of  H.  emissions  in  the  CH4 
discharges  with  that  of  the  much  better  characterized  H2  discharges1’ 2' 14-16  provides  the 
desired  H-atom  production  estimate.  The  correlation  of  bright  C2  band  emission  with 
discharge  inhomogeneity  is  also  documented  and  discussed.  The  issue  of  discharge 
homogeneity  may  be  important  primarily  from  the  viewpoint  of  further  studies  involving 
laser  spectroscopic  measurements  of  radical  fragment  yield  as  a  function  of  E/N,  and 
determination  of  scaling  parameters. 


Experimental  Details 

A  schematic  diagram  of  the  experimental  apparatus  is  shown  in  Figure  15.  The 
main  components  of  the  apparatus  are  the  discharge  cell,  the  light  detection  electronics, 
and  the  high  voltage  pulser.  The  discharge  cell  was  a  cylindrical  glass  enclosure, 
approximately  37  cm  in  height  and  7.5  cm  in  diameter.  Two  sets  of  opposed  fused  silica 
windows  (5-cm  diameter),  orthogonal  to  each  other,  were  available  as  viewports.  The 
cell  was  evacuated  through  3. 5-cm  i.d.  glass  and  metal  lines  by  a  roughing  pump.  The 
base  pressure  for  the  cell  was  3  mTorr. 

The  discharge  gas  flow  was  metered  by  flow  controllers  (MKS  Type  1259  flow 
controller  and  a  MKS  Type  247C  four-channel  readout),  while  the  cell-pressure  was 
governed  by  a  control  valve  and  an  exhaust  valve  controller  (MKS  Types  253A  and  252, 
respectively).  The  gas  pressure  was  measured  with  a  baratron  pressure  gauge  (MKS 
Type  102  A)  whose  output  was  sent  to  the  valve  controller  to  maintain  constant  pressure 
in  the  discharge  cell.  The  gas  flow  could  be  varied  between  0  and  100  SCCM,  and  the 
baratron  operated  in  the  pressure  range  0-100  Torr.  Gas  pressure  and  flow  rate  were  set 
by  input  to  the  PC.  Throughout  this  study  the  gas  flow  rate  was  fixed  at  30  SCCM  and 
the  pressure  (p)  at  50  Torr.  Discharge  gases  were  UHP  grade  (99.999%). 

The  discharge  electrodes  were  two  aluminum  discs,  3.2  cm  in  diameter  with  an 
interelectrode  spacing  (d)  that  varied  between  0.25  -  1.0  cm  in  this  work  (pd  range  12.5  - 
50  Torr  cm).  The  discharge  pulse  repetition  rate  used  was  10  Hz  throughout.  The 
interelectrode  space  could  be  illuminated  by  an  UV  preioniaztion  lamp.  The  high  voltage 
pulser  constructed  for  this  study  powered  the  discharge  and  is  based  upon  a  stacked 
MOSFET  switch.  Output  current/voltage  measurements  were  made  using  high  voltage 
probes  (HP  Model  P6015)  across  the  discharge  and  across  a  series  current  sense  resistor. 
The  probe  outputs  were  monitored  by  a  digital  oscilloscope  (EeCroy  9354C)  for  signal 
capture  and  averaging  (100  shots).  Current  and  voltage  waveforms  were  then  transferred 
to  a  PC  for  processing. 

High  speed  direct  discharge  imaging  and  spectral  imaging  were  done  using  a 
gateable,  intensified  CCD  camera  (Princeton  Instruments  Model  ICCD-576).  Spectral 
images  were  made  using  10-nm  bandpass  filters  at  43 1  nm  (CH)  and  516  nm  (C2).  Time 
resolved  H  ,  CH  or  C2  emissions  from  the  discharge  were  detected  using  a  gated 
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photomultiplier  tube  (PMT,  Hamamatsu  R1477).  Light  was  collected  through  an  iris 
diaphragm  by  a  focusing  lens  and  passed  through  the  various  bandpass  filters  onto  the 
PMT.  The  PMT  was  under  a  constant  bias  throughout  the  experiment,  and  was  gated  on 
for  a  5-  s  period  encompassing  the  discharge  pulse  duration.  Neutral  density  filters  were 
used  to  avoid  saturation.  The  PMT  output  was  amplified  and  sent  to  a  digital 
oscilloscope  for  signal  capture  and  averaging  (100  shots),  and  to  a  PC  as  above.  Single¬ 
shot  C2  emission  profiles  were  also  recorded.  When  discharge  imaging  is  not  being 
conducted,  a  second  gated  PMT  in  place  of  the  camera  allows  for  simultaneous  time 
resolved  emission  spectroscopy  of  a  pair  of  methane  fragments. 

Results 

Discharge  Characterization 

Averaged  current,  gap  voltage  and  power  waveforms  for  discharges  in  pure 
methane  (flows)  at  50  Torr  and  using  a  5  kV  pulser  drive  voltage  are  shown  in  Figure  16 
for  four  different  electrode  spacings  (1.0,  0.75,  0.50  and  0.25  cm).  The  pulser  provides 
rapid  excitation  pulse  rise  and  fall  times  (typically  <  20  ns),  control  over  pulse  duration, 
and  low  jitter  (typically  <10  ns).  Note  that  the  discharge  current  turns  on  earlier  as  the 
gap  is  decreased.  A  fixed  trigger  pulse  width  of  440  ns  was  used  throughout  this  work, 
(and  the  discharge  pulse  width  variation  results  from  that  turn-on  trend.) 

The  current  pulse  rise  is  typically  delayed  relative  to  that  of  the  gap  voltage  pulse 
by  roughly  25  ns  in  each  case,  and  there  is  a  substantial  overlap  between  current  and 
voltage  during  the  short  initial  peak  in  the  voltage  pulse.  The  power  curves  reflect  the 
degree  of  overlap  with  the  brief  (~30  ns  FWHM),  high-power  initial  peak  that  occurs  at 
each  gap.  Additionally,  no  indications  of  arc  formation  are  seen  by  way  of  sudden 
voltage  drops  and  associated  current  excursions;  pulse-to-pulse  shape  variations  are 
minimal. 

In  all  cases  the  gap  voltage  falls  to  a  lower  steady  state  value  after  the  initial  peak 
until  the  pulser  is  triggered  off.  This  implies  that  discharge  E/N  is  substantially  higher 
during  the  initial  pulse  peak  than  in  the  relatively  long  final  steady  state  period.  In  order 
to  determine  E/N  from  the  gap  voltage  and  the  discharge  conditions  (50  Torr,  300  K),  the 
voltage  drop  across  the  thin  (--50  m)  cathode  sheath  must  subtracted.  The  sheath  drop 
and  average  E/N  for  the  discharges  were  estimated  during  the  steady  state  portion  by 
plotting  the  average  voltage  (between  200  -  400  ns)  against  electrode  separation,  as 
shown  in  Figure  17.  Though  pd  varies  by  a  factor  of  4,  the  data  bear  the  usual  linear 
relationship  ’  ,  and  indicates  reasonable  values  for  both  the  sheath  drop  of  300  V,  and 
an  average  E/N  of  80  Td.  When  that  sheath  voltage  drop  is  subtracted  from  the  measured 
gap  voltages,  shown  in  Figure  16,  average  E/N  during  the  steady  state  portion  of  the 
discharges  agrees  well  with  that  value.  E/N  conditions  during  the  initial  breakdown 
portion  of  the  discharge  pulses  range  up  to  a  factor  of  four  higher  than  the  steady  state 
portion  of  the  pulse. 
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The  range  of  E/N  values  found  here  is  attractive  from  the  viewpoint  of  assessing 
the  hydrocarbon  fragment  yield  as  a  function  of  E/N.  To  determine  if  the  discharge  can 
also  be  made  sufficiently  diffuse  to  facilitate  fragment  yield  measurements  by  laser 
spectroscopy  techniques,  discharge  imaging,  spectral  imaging,  and  time-resolved 
spectroscopy  techniques  were  used.  Due  to  the  short  duration  of  a  typical  discharge 
pulse,  the  discharges  were  imaged  with  a  gateable,  intensified  CCD  camera  to  determine 
the  degree  of  homogeneity  and  pulse-to-pulse  variability. 

Discharge  V-I  pulse  shapes,  especially  risetimes,  and  discharge  pd  values,  are 
important  detenninants  of  discharge  uniformity17’ 18,  as  are  levels  and  methods  of 
discharge  gas  preionization19.  Using  discharge  conditions  identical  to  those  shown  in 
Figures  16  and  17,  Figure  18  shows  typical  100  ns  exposures  documenting  the  (dramatic) 
effect  that  decreasing  the  electrode  spacing  has  upon  discharge  uniformity.  The 
transition  between  multiple  streamer  discharge  and  diffuse  discharge  occurs  under  the 
indicated  conditions  between  0.5  cm  and  0.25  cm  gaps,  or  below  a  pressure-discharge 
gap  product  (pd)  of  25  Torr-cm.  (All  discharge  gaps  studied  show  sufficient 
homogeneity  to  allow  for  the  application  of  laser  spectroscopy  techniques.) 

The  behavior  of  the  discharge  as  a  function  of  pulser  drive  voltage  is  shown  in 
Figure  19  for  a  fixed  gap  of  1.0  cm  and  drive  voltages  of  4,  6  and  8  kV.  All  other 
experimental  parameters  are  identical  to  those  discussed  above.  Features  here  are  similar 
to  those  of  Figure  16,  except  that  the  current  risetime  and  turn  on  time  are  much  slower 
for  the  lowest  drive  voltage.  Here  too  the  gap  voltage  pulses  peak  just  before  the 
corresponding  current  turns  on,  the  degree  of  overlap  being  reflected  in  the  corresponding 
power  curves.  The  overvoltage  at  the  beginning  of  the  pulse  is  clearly  most  effective  for 
power  deposition  at  high  E/N  conditions. 

Extracting  the  temporal  evolution  of  discharge  E/N  from  the  above  gap  voltage 
pulses  was  accomplished  by  subtracting  a  scaled  sheath  voltage  drop  from  the  gap 
voltage  and  accounting  for  discharge  conditions  (300  K,  50  Torr,  1  cm  gap).  The  result  is 
shown  in  Figure  20.  Again  note  the  overlap  with  the  corresponding  current  pulses:  the 
high  E/N  values  early-on  for  the  4  kV  drive  voltage  is  largely  ineffective  given  the  late 
turn  on  time  for  that  case.  However,  the  high  peak  power  for  the  higher  drive  voltages 
substantially  overlaps  the  early  high  E/N  values,  a  result  that  is  a  major  determinant  of 
time  resolved  hydrocarbon  fragment  emission  profiles. 

Time  resolved  H  and  CH  band  emissions  (43 1  mn)  for  experimental  conditions 
identical  to  those  shown  in  the  previous  two  figures  are  shown  in  Figure  2 1 .  The  effect 
on  direct  electron  impact  excited  methane  fragment  (H(n=3)  and  CH(A2  .  states) 
production  of  the  early  high  E/N  values  for  the  higher  drive  voltages  is  clear,  though  the 
emission  profiles  and  relative  intensities  are  somewhat  different.  Both  emissions  rise  at 
the  same  rate  as  the  discharge  power  pulse  indicating  direct  electron  impact  dissociative 
excitation  of  methane  as  their  sources.  The  brighter  (factor  of  five)  CH  emission  has  a 
markedly  long  fall  time  due  to  the  long  radiative  lifetime  and  the  low  quenching"  rate 
compared  to  very  strongly  quenched"  ’  H(n=3)  atomic  level.  At  50  Torr  the  H 

emission  profile  very  closely  resembles  that  of  the  discharge  power  for  all  drive  voltages. 
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For  the  higher  drive  voltages  there  is  a  short  current  peak  at  the  beginning  of  the  pulses. 
However,  both  the  H  and  CH  emissions  are  much  more  significantly  enhanced  at  the 
beginning  of  the  pulse  than  current  pulse  shapes  would  suggest,  thus  showing  the  effect 
on  methane  fragmentation  which  results  from  the  high  E/N  conditions  at  the  early  on. 

The  final  time  resolved  emission  results  from  CH4  discharges  reported  here  are 
for  the  C2  Swan  bands  near  516  mn.  Complemented  by  (direct)  discharge  imaging 
and  spectral  imaging  of  the  C2  and  CH  bands,  the  rare  occurrence  of  highly 
inhomogeneous  discharges  can  be  partially  understood.  On  rare  occasions  (direct) 
discharge  images  indicated  the  presence  of  only  one  or  two  very  bright  streamers  in  the 
discharge.  Almost  entirely  absent  at  the  smallest  gap  studied,  the  incidence  of  such 
highly  inhomogeneous  discharge  pulses  grew  to  roughly  3  percent  of  pulses  at  the  largest 
gap  of  1.0  cm.  Further,  a  dramatic  increase  in  bright  streamer  formation  frequency 
occurred  in  the  absence  of  irradiation  of  the  discharge  space  by  the  UV  preionization 
source.  Insights  into  a  potential  cause  of  the  formation  of  a  few  bright  streamers  as 
opposed  to  multiple  streamer  discharges  were  suggested  by  related  on-going 
investigations  .  Vibrationally  hot  C2  bands  have  been  observed  in  similar  discharges, 
and  at  higher  pressures  (100  Torr)  visible  carbon  deposits  have  been  observed  to  fonn  on 
the  discharge  electrodes. 

•5 

Figure  22  shows  five  single-shot  current  pulses  and  their  corresponding  C2  (A  g 
-  X3  u,  v  =  0)  emission  profiles  (near  516  mn),  taken  with  and  without  UV 
preionization  of  the  discharge  (gas).  (Direct)  Images  of  those  same  pulses  are  shown  in 
Figure  23.  The  discharge  conditions  were  0.75  cm  discharge  gap,  50  Torr  methane,  and 
drive  voltage  of  5  kV.  Except  for  a  delayed  single  shot  in  the  case  without  UV 
irradiation,  the  current  pulses  are  essentially  identical  between  the  two  cases,  and  clearly 
quite  repeatable.  The  clearest  difference  between  the  two  data  sets  is  the  much  greater 
brightness  without  UV  irradiation.  Inspection  of  Figure  23  correlates  the  bright  streamer 
discharge  conditions  with  bright  C2  emission.  The  rise  of  the  more  typical  diffuse,  weak 
emissions  in  the  case  with  UV  irradiation  essentially  follows  the  current.  (This  suggests 
excitation  by  direct  electron  impact,  requiring  the  carbon  molecules  to  be  normal,  low 
density  gap  constituents.) 

Clearly,  the  discharge  fully  dissociates  some  fraction  of  methane,  and  some 
carbon  atoms  subsequently  combine  to  form  C2  and  presumably  larger  carbon  clusters. 
The  rise  in  the  bright  carbon  emissions  associated  with  bright  streamer  appearance  is 
longer  and  slower  to  steady  state,  after  an  initial  rapid  rise  similar  to  the  weak  emissions 
case.  The  brighter  emission  and  different  temporal  emission  profile  suggest  the  possible 
involvement  of  large  carbon  clusters  breaking  up  under  direct  electron  impact  within  the 
gap,  or  being  detached  from  the  cathode  and  breaking  up  as  the  pulse  rises.  Carbon 
clusters  that  formed  during  repetitive  pulsing  of  the  discharge  could 
adhere  to  the  electrodes  providing  both  the  source  of  the  bright  streamers  as  well  as  the 
vibrationally  hot  C2  bands.23  Such  processes  would  locally  enhance  the  discharge 
conductivity  and  contribute  to  the  fonnation  of  inhomogeneous  plasmas  with  one  or  two 
bright  streamers.  The  C2  emission  profiles  in  the  case  of  highly  inhomogeneous  bright 
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streamer  discharges  are  different  from  any  other  observed  emission  profile  (see  Figure 
21),  suggesting  different  kinetics  involved  in  excitation  of  the  C2  molecule  in  those  cases. 

The  final  experimental  results  associated  with  the  discharge  characterization  show 
the  spatial  distribution  of  C2  emissions  in  comparison  with  that  of  the  brightest 
hydrocarbon  fragment  emission  bands  of  CH.  Figure  24  shows  images  of  the  discharges 
similar  to  those  of  Figure  23,  except  that  they  have  been  imaged  through  bandpass  (10 
nm  FWHM)  filters.  The  typical  spatial  distribution  of  the  C2  emission  is  quite  uniform 
and  weak  (upper  left)  when  the  discharge  is  diffuse,  indicating  a  largely  uniform 
distribution  of  C2.  Bright  streamer  formation  is  clearly  accompanied  by  strongly 
enhanced  C2  emission  (upper  right).  The  spatial  distribution  of  the  CH  emission  is  quite 
similar  to  that  of  C2  regardless  of  the  degree  of  discharge  homogeneity. 


H-atom  Production 


Using  the  technique  of  actinometry24'27,  an  estimate  of  H-atom  production  during 
a  discharge  pulse  can  be  made  by  comparing  the  H  emission  from  the  methane 
discharge  (Figure  21)  to  that  from  the  very  well  characterized  pure  H2  discharge1’2' 14-16 
under  the  same  experimental  conditions.  Gap  voltage,  discharge  power,  and  the  resulting 
H  emission  from  discharges  into  pure  hydrogen  at  50  Torr  using  a  6  k V  pulser  drive 
voltage  are  shown  in  Figure  25.  The  temporal  behavior  of  the  gap  voltage  and  discharge 
power  are  very  similar  to  the  analogous  case  for  discharges  into  methane  (see  Figure  19, 

6  kV  drive  voltage  case),  although  both  fall  more  rapidly  in  the  hydrogen  plasma  during 
the  latter  portion  of  the  discharge  pulse.  This  behavior  clearly  indicates  that  discharge 
E/N  in  the  pure  hydrogen  discharges  decreases  somewhat  more  during  the  latter  portion 
of  the  discharge  pulse,  as  opposed  to  remaining  essentially  constant  after  the  initial  high 
E/N  peak  in  methane  (see  Figure  20,  6  kV  case).  The  likely  cause  of  the  minor  difference 
in  temporal  behavior  is  that  the  methane  plasmas  are  slightly  electronegative.  The  H 
emission  from  the  pure  hydrogen  discharge  shown  in  Figure  25  follows  the  discharge 
power  pulse,  identical  to  the  temporal  behavior  of  the  H  emission  from  methane 
(compare  to  Figure  21,6  kV  case). 

Discharges  in  pure  hydrogen  have  been  extensively  studied1'2'14-16,  and  the 
fraction  of  the  discharge  power  deposited  into  the  various  inelastic  modes  of  hydrogen  as 
a  function  of  E/N  has  been  established.  For  E/N  above  80  Td  the  fraction  of  the 
discharge  power  producing  dissociation"  of  H2  into  two  H  atoms  is  greater  than  50  %. 
Absolute  H  atom  density  and  H  emission  measurements1  for  such  discharges  have  also 
made  as  a  function  of  discharge  E/N.  Thus,  it  is  possible  to  write  the  number  density  of 
H  atoms  in  pure  hydrogen  discharges  as, 


2^X(E/N)  PH2dt 


D  V 


(1) 


where  Pm  is  discharge  power,  X(E/N)  is  the  known  fractional  power  deposition  into 
dissociation2  as  a  function  of  E/N,  De  is  the  dissociation  energy  of  H2  (9.8  eV),  and  V  is 
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the  discharge  volume.  The  factor  of  2  accounts  for  the  fact  that  each  dissociation  event 
produces  two  hydrogen  atoms.  Equation  (1)  can  be  used  to  directly  determine  the  H- 
atom  production  in  pure  hydrogen  discharges  from  electrical  measurements  alone. 

H  emission  from  such  a  discharge  results  from  the  competition  between 
production  of  H(n=3)  atoms,  and  losses  by  radiation,  or  by  collisional  quenching,  to  the 
ground  state.  As  noted  above,  the  rapid  rise  of  H  emission  (in  concert  with  the  rise  in 
the  discharge  power  pulse)  indicates  that  production  of  H(n=3)  atoms  in  this  study 
proceeds  by  direct  electron  impact  dissociative  excitation  of  the  parent  gas.  The  resulting 
H  emission  intensity  can  be  written”  '  as  the  ratio  of  the  rate  of  production  to  that  of 
loss: 


H2  Chv,2A,2k»2<2\H2\ 

Ha~  X^-Cl//,l  '  U 

j 

Here  nf 2  is  the  electron  density  in  the  hydrogen  discharge,  h  32  and  A32  are  the  photon 
energy  and  Einstein  coefficient  for  spontaneous  emission  for  the  transition,  respectively, 
Ay  is  the  sum  of  the  A  coefficients  for  all  radiative  de-excitation  processes  from  the 
H(n=3)  level  (H  and  Lyman  ),  and  C  is  a  constant  that  contains  all  optical  and 
geometrical  parameters  such  as  solid  angle  viewed,  transmissions,  and  detection  system 
response.  The  coefficients  kf2 and  kg2  quantify  electron  impact  dissociative  excitation 

of  H2  to  H(n=3)  +  H(ls)  and  the  quenching  rate  of  H(n=3)  atoms  in  hydrogen, 
respectively.  The  square  brackets  denote  concentration. 


For  discharges  into  pure  CH4,  the  resulting  H 
written  as, 


JX  +  tJTtcwj 

j 


emission  intensity  can  similarly  be 


(3) 


where  the  symbols  have  analogous  meanings.  The  ratio  of  H 
discharges  then  is  simply, 


jH 2 
1  Ha 
tCHA 
Ha 


J 


emissions  from  the  two 


(4) 


It  has  been  shown  that  H-atom  production  in  (fast)  pulsed  discharges  in  hydrogen 
(measured  by  two  photon  allowed  laser  induced  fluorescence  or  TALIF)  is  directly 
proportional  to  the  H  emission  signal1 .  The  same  relation  between  H-atom  production 
and  H  emission  should  hold  for  discharges  in  methane  since  the  thresholds  for  both 
electron  impact  dissociative  excitation  ’  (DE)  and  dissociation^  ’  (D)  are  nearly  the 
same  in  both  hydrogen  and  methane,  i.e. 


(E“2(DE)  =  19  eV,  E™\DE)  =  20  eV,  and  E“2  (D)  =  9.8 eV,E™\D)  =  8.8eV). 


The  threshold  for  direct  electron  impact  excitation  of  H  atoms  to  H(n=3)  is,  of 
course,  independent  of  the  parent  gas.  Together,  these  facts  allow  replacement  of  the 
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ratio  of  production  rates  of  H(n=3)  atoms  by  dissociative  excitation  in  Equation  (4)  with 
the  corresponding  rates  of  direct  electron  impact  excitation  of  ground  state  H  atoms  to  the 
H(n=3)  level.  With  that  substitution  Equation  4  becomes, 


/ 


H  2 
Ha 


or™ 


HI 


tCH  4 
1  Ha 


„  CH  4  CH  4 


[H] 


CH  4 


(5) 


where  the  subscript  ei  indicates  the  rate  coefficient  for  electron  impact  of  H  atoms  to 
H(n=3),  and  [II] m  and  [H]ch4  refer  to  the  concentration  of  H  atoms  in  the  hydrogen  and 
methane  discharges,  respectively. 


The  product  of  the  direct  electron  impact  excitation  rate  coefficient  and  electron 
density  in  the  each  plasma  is  proportional  to  the  fraction  of  discharge  power  devoted  to 
that  process  multiplied  by  the  power  deposited  in  the  gas,  i.e. : 


H  2  H  2 


oc  Y(E  !  N)  PH 2 ,  and  k 


CH4  CH  4 


oc  Y'(E  /  N)  Pc 


CH  4’ 


(6) 


2  • 

where  Y(E  /  N)  is  the  known"  fractional  power  deposited  into  electron  impact  production 
of  H(n=3)  as  a  function  of  E/N  in  hydrogen,  and  Y  (E /  N)  is  the  analogous  factor  in 


methane  discharges.  Note  that  the  proportionality  constants  in  Equations  (6)  are  the 
inverse  of  the  excitation  potential  of  H(n=3),  and  are  therefore  identical  for  each  of  those 
relations.  Substitution  of  Equations  ( 1 )  and  (6)  into  Equation  (5),  and  assuming  Y  «  7' , 
the  final  expression  used  to  estimate  H  atom  production  in  the  CH4  discharges  of  interest 
can  be  written  as, 


[H] 


T  CH  A  p 
1  Ha  rH  2 


CH  4 


[H2 
!  Ha 


CH  4 


YjAsi+ka  [Ctf-il  2 ^X(EIN)PH2dt 

YaPPiP]  Ev 
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(V) 


In  this  relation,  the  relevant  quenching  rate  constants  for  H(n=3)  atoms2122  are  2  x  10'9 
cm  /s  in  H2  and  3.5  x  10'  cm/s  in  CH4,  and  the  measured  H  intensity  profiles  needed  as 
inputs  are  shown  in  Figures  2 1  (6  kV  drive  voltage  case)  and  25  for  emission  from 
methane  and  hydrogen  discharges,  respectively.  The  discharge  power  profiles  are  shown 
in  Figure  19  for  methane  (6  kV  drive  voltage  case)  and  Figure  25  for  hydrogen.  The 
factor  X was  taken  to  be  0.6  during  the  initial  high  E/N  peak  in  the  power  pulse,  and  0.45 
during  the  much  longer  steady  state  portion  of  the  discharge  (see  Figures  20  and  25). 


Using  Equation  (1),  the  H-atom  production  during  a  single  discharge  pulse  into  50 
Torr  of  pure  hydrogen  using  a  pulser  drive  voltage  of  6  kV  was  calculated  and  is  shown 
as  the  upper  plot  in  Figure  26.  Similarly,  using  Equation  (7),  the  H-atom  production 
during  a  single  discharge  pulse  into  methane  under  the  same  experimental  conditions  is 
shown  in  the  lower  plot  of  that  figure.  For  comparison  purposes,  the  discharge  pulse 
power  and  the  associated  emission  are  also  shown  in  that  plot.  Note  that  the 
production  rate  is  significantly  higher  during  the  brief  high  E/N  early  peak  in  the 
discharge  pulse  that  is  typical  in  such  fast  discharges.  Production  of  =3  x  10 14  H  atoms 
per  cubic  centimeter  of  discharge  volume  during  a  400  ns  discharge  pulse  is  significant. 
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The  total  energy  deposited  in  the  gas  during  a  pulse  is  5.6  mJ,  and  the  energy  deposited 
per  unit  volume  is  0.7  mJ/cm  .  Verification  of  that  level  of  H-atom  production  will 
require  the  use  of  (more  elaborate)  two-photon  laser  induced  fluorescence  or  mass 
spectrometric  measurements  CH4  fragments. 


Discussion  and  Conclusions 

Techniques  for  the  control  of  the  transition  from  diffuse  to  inhomogeneous 
plasmas  in  high  E/N  hydrocarbon  discharges  have  been  established.  (Plasma 
homogeneity  can  be  made  sufficient  good  to  allow  for  further  quantitative  study  of  the 
hydrocarbon  fragment  yield  by  laser  spectroscopic  techniques).  The  steady  state  portion 
of  the  discharge  voltage  pulse  follows  the  classic  pd  scaling17'18  relationship  over  the  pd 
range  studied  (12.5  -  50  Torr  cm).  The  resulting  plasmas  are  typically  diffuse  in  the  lower 
end  of  that  pd  range,  and  even  at  the  higher  pd  values  are  sufficiently  uniform  (multiple 
streamers),  (for  laser  diagnostic  techniques  to  be  applied.) 

The  association  of  bright  C?  emission  with  the  most  inhomogeneous  discharge 
plasmas,  and  the  strong  inhibiting  influence  that  discharge  gas  preionization  has  over 
their  formation,  are  clearly  established.  The  possible  role  of  carbon  clusters  in  the 
formation  of  plasma  inhomogeneities  should  be  studied  further.  A  potentially  attractive 
method  for  carbon  concentration  modification  in  the  discharges  could  be  the  addition  of 
N2.  Any  eventual  scramjet  ignition  system  would  have  to  allow  for  the  presence  of 
atmospheric  nitrogen,  and  the  fonnation  of  CN  radicals  should  result  in  a  reduction  in  the 
carbon  concentration  and  attendant  cluster  formation.  CN  formation  (should)  can  be 
easily  monitored  by  optical  spectroscopy. 

Significant  levels  of  hydrogen  atom  production  by  fast,  high  E/N  discharges  are 
estimated  from  the  measured  H  emission  profiles  and  by  the  application  of  an 
actinometric  model.  The  results  are  quite  encouraging  from  the  viewpoint  of  application 
to  scramjet  ignition  and  combustion  sustainment.  (Verification  of  that  estimate  will 
require  much  more  elaborate  experimental  investigations.) 
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Figure  15.  Schematic  diagram  of  the  experimental  apparatus. 
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Gap  Voltage  (kV) 


Figure  16.  Current,  gap  voltage,  and  power  characteristics  for  discharges  into  50 
Torr  CH4  and  5  kV  pulser  drive  voltage  as  functions  of  discharge  gap. 
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Electrode  Separation  (cm) 


Figure  17.  Average  gap  voltage  between  200  -  400  ns  plotted  against  anode-cathode  gap 
(see  Figure  16). 
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5  kV,  50  Torr,  0.5  cm  gap 


Is 

5  kV,  50  Torr,  0.25  cm  gap 


Figure  18.  Direct,  high-speed  two-dimensional  CH4  discharge  images  as  a  function 

of  interelectrode  spacing. 
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Figure  19.  Current,  gap  voltage,  and  power  characteristics  as  functions  of  pulser 
drive  voltage  for  discharges  into  50  Torr  CH4  and  a  1.0  cm  electrode  gap. 
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Figure  20.  Temporal  evolution  of  the  discharge  E/N  for  three  drive  voltages,  a  1  cm 
gap,  and  50  Torr  methane. 
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Relative  CH  Intensity  Relative  H-alpha  Intensity 


Time  (ns) 


Figure  2 1 .  Temporal  evolution  of  the  most  intense  visible  emissions  from  the  excited 
hydrocarbon  fragments:  H  and  CH  (43 1  mn  band).  Note  that  the  H  emission 
essentially  follows  the  discharge  power  pulse  while  the  CH  emission  peak  occurs  later 
in  time.  Experimental  conditions  are  the  same  as  for  the  previous  two  figures. 
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Relative  Intensity 


Time  (ns) 


Time  (ns) 


Figure  22.  Single  shot  current  pulses  and  the  corresponding  C2  Swan  band  emission  profdes 
near  516  nm  for  cases  with  and  without  UV  irradiation  of  the  discharge  gas.  Discharge 
conditions  are  50  Torr  methane,  5  kV  drive  voltage  and  fixed  gap  of  0.75  cm.  The 
corresponding  direct  discharge  images  are  presented  in  Figure  23. 
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Figure  23.  Direct  discharge  images  of  the  same  discharge  pulses  shown  in  Figure  24,  with 
UV-irradiated  discharges  on  the  left.  Bright  streamer  formation  is  much  more  common 

when  the  discharge  gas  is  not  preionized. 
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typical  carbon  dimer  emission 


intense  carbon  dimer  emission 


(with  UV  preionization) 


(no  UV  preionization) 


CH  radical  emission 


(no  UV  preionization) 


Figure  24.  Spatial  distribution  of  C2  and  CH  emissions  for  typical  discharges  (left),  and 
highly  inhomogeneous  bright  streamer  discharges  (right). 
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Figure  25.  Current  and  gap  voltage  characteristics  (top)  and  the  resulting  Ha  emission 
(bottom)  for  discharges  into  pure  H2  at  50  Torr  pressure,  1.0  cm  gap,  and  a  drive  voltage  of  6 
kV  (each  area  100  shot  averages.  Compare  to  Figures  19  and  21  for  CH4  discharges. 
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Figure  26.  Estimated  H-atom  production  in  pure  hydrogen  (top)  and  methane 
(bottom),  each  during  a  discharge  pulse.  Both  are  determined  from  experiments  using 
a  1  cm  gap,  a  pulser  drive  voltage  of  6  kY,  and  50  Torr  gas  pressure. 
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